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Abstract 

In the study of chaotic behaviour of systems of many hard spheres, Lyapunov 
exponents of small absolute value exhibit interesting characteristics leading to spec- 
ulations about connections to non-equilibrium statistical mechanics. Analytical ap- 
proaches to these exponents so far can be divided into two groups, macroscopically 
oriented approaches, using kinetic theory or hydrodynamics, and more microscopi- 
cally oriented random-matrix approaches in quasi-one-dimensional systems. In this 
paper, I present an approach using random matrices and weak disorder expansion in 
an arbitrary number of dimensions. Correlations between subsequent collisions of a 
particle are taken into account. It is shown that the results are identical to those of a 
previous approach based on an extended Enskog-equation. I conclude that each ap- 
proach has its merits, and provides different insights into the approximations made, 
which include the Stofizahlansatz, the continuum limit, and the long- wavelength 
approximation. The comparison also gives insight into possible connections between 
Lyapunov exponents and fluctuations. 
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1 Introduction 



In recent years, investigations into the connections between the theory of dy- 
namical systems and non-equihbrium statistical mechanics have yielded many 
interesting and important results. Gallavotti and Cohen [US], for instance, 
conjectured that many-particle systems as studied by statistical mechanics 
will generally be strongly chaotic. This has prompted a great deal of interest 
in the connections between chaos on the one hand and the decay to equilibrium 
and transport coefficients on the other (see for instance Ref. [3]). A central 
role in the study of chaos and related properties is played by the Lyapunov 
exponents, which describe the exponential divergence or convergence of nearby 
trajectories in phase space. 

Some of this interest has been directed towards the Lyapunov exponents of 
the propotype system of many hard spheres. Several analytical calculations 
of, among other things, the largest Lyapunov exponent and the sum of all 
positive Lyapunov exponents have been performed [H[5]|6|7|8|9] . Lyapunov 
exponents of many-particle systems have also been evaluated numerically in 
molecular-dynamics simulations (see, for instance Ref. [TO|ITT] ). Because of 
their unexpected behaviour, in particular the Lyapunov exponents of small 
but nonzero absolute value have received attention. A step structure occurs 
in the Lyapunov spectrum near zero whenever the system is large enough 
compared to the mean free path, as was first noted by Posch and Hirschl [12] 
and later also found in other systems (see, for example, Refs. p!T|T3] ). These 
Lyapunov exponents differ from the exponents of larger absolute value, in the 
sense that all particles contribute to them, much like in the case of the zero 
Lyapunov exponents, and the corresponding modes appear to be, on average 
and to first approximation, linear combinations of these zero modes with a 
sinusoidal modulation in the position. Initially, it was hoped that describing 
the Lyapunov modes through a macroscopically oriented approach such as 
hydrodynamics or an Enskog equation might provide insight into possible 
connections between chaos and transport. In Ref. [7], it has been shown that 
the small exponents can in fact be viewed to belong to Goldstone modes and 
that the behaviour found in simulations [12] can be understood from this. A set 
of equations was derived for these exponents by the use of an extended Enskog 
equation and values for the exponents were obtained. Other attempts to 
understand these exponents have been based on hydrodynamic equations [14J , 
and, although limited to quasi-one- dimensional systems, random matrices 
along with approximations of weak disorder [T5|16|17j . 

In view of the two distinct approaches to the Goldstone modes, through ran- 
dom matrices on the one hand and through the Enskog equation on the other. 



and his continued interest. 
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it is of interest to investigate whether the results of Ref. [7] can also be derived 
using techniques from random matrix theory. In this paper, instead of starting 
from the Enskog equation, I make use of random matrices and the weak- 
disorder expansion. Unlike the previous random-matrix approaches mentioned 
above, the present derivation is not limited to quasi-one-dimensional systems. 
The approximations needed to arrive at quantitative results can be studied 
more carefully in some cases, and are similar to those used in the derivation of 
Ref. [7j. By comparing the Enskog and random-matrix approaches, one gains 
insight into the approximations made in both approaches and the associated 
inaccuracies. Of special interest are the consequences of the thermodynamic 
limit, since finite-size effects in the Lyapunov exponents may be related to 
fluctuations and decay of correlations. 

This paper is organised as follows. In Sec. |2l Lyapunov exponents are briefly 
introduced as well as the dynamics in tangent space of freely moving hard 
spheres in tangent space. Next, in Sec. [3l a summary is given of the Goldstone 
modes and the calculation of Ref. [7] by the use of an extended Enskog 
equation. In Sees. IH El and O it is explained how the results found from 
the extended Enskog equation can also be derived through the use of random 
matrices. The approaches are compared in Sec. [TJ and the approximations 
needed are discussed. Possible corrections are considered and it is pointed out 
how these may lead to insight in the connections between non-equilibrium 
behaviour and chaotic properties. 



2 Lyapunov exponents and the dynamics in tangent space 

Consider a d-dimensional system of particles moving in a 2(iA^-dimensional 
phase space F. At time t = 0, the system is assumed to be in an initial point 
7o in this phase space, from which it evolves with time according to 7(70, t). If 
the initial conditions are perturbed infinitesimally by ^70, the system evolves 
along an infinitesimally different path 7(70, t) -|- 57(70,^), where 6'j denotes a 
coordinate in the tangent space 6T, and 57(70, 0) = ^70. The evolution of a 
vector in the tangent space is described by 



57(70, t) = M^,(i) ■'^70 , 




where M^(,(t) is a 2 dA^- dimensional matrix defined by 



(2) 
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The Lyapunov exponents are the possible average asymptotic growth rates 
of infinitesimal perturbations 57^(7, t) associated with the eigenvalues of 
M^„(t), i.e., 



\i = lim - (In I /ij (t) I + z arg /Xj (t ) ) 



(3) 



If the system is ergodic, it will eventually come arbitrarily close to any point 
in phase space for all initial conditions except for a set of measure zero. The 
Lyapunov exponents are thus the same for almost all initial conditions. In the 
literature, one also finds the Lyapunov exponents defined with reference to 
the eigenvalues of [M-yo(t)^ ■ M^Q(t)]2, in which case they are real. 

The symmetries of the dynamics of the system generate vectors in tangent 
space which do not grow or shrink exponentially and therefore have Lyapunov 
exponents equal to zero. For a system of hard spheres under periodic boundary 
conditions, these symmetries and their corresponding zero modes are uniform 
translations, Galilei transformations, time translations, and velocity scaling. 

We now consider a gas of identical hard spheres of diameter a and mass m 
in d dimensions in the absense of external fields. As there are no internal 
degrees of freedom, the phase space may be represented by the positions and 
velocities Vj of all particles, enumerated by i, and similarly the tangent space 
by infinitesimal deviations 5rj and 5vj. The evolution of the system in phase 
space consists of a sequence of free fiights interrupted by collisions. During the 
free fiights, the particles do not interact and their positions change linearly 
with the velocities; similarly, 6r changes linearly with 5v. For rigid spheres 
the collisions are instantaneous. At the moment of the collision, momentum 
is exchanged between the two particles involved along the collision normal 
CT = (rj — rj)/a at impact, as shown in Fig. [H At the instant of the collision, 
none of the other particles are assumed to interact. 

From Eq. ([2]) and the phase-space dynamics, the dynamics in tangent space 
can be derived [HH]. During the free fiight between the instant of a collision 
tz {z being the number of the collision in the sequence) and t, there is no 
interaction between the particles and the components of the tangent-space 
vector transform according to 




(4) 



(5) 
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Fig. 1. Geometry of a collision of two particles i and j of diameter a, in relative 
position Tij = Tj — r,- and with relative velocity v- 



j — Vj. The collision normal 
<T is the unit vector pointing from the centre of particle j to the centre of particle 
i. The circle drawn represents the locus of closest approach. Coordinates before 
the collision are marked with a prime. The dashed lines indicate an infinitesimally 
displaced path. 

in which I is the d x d identity matrix. 

As shown in Fig. [U infinitesimal differences in the positions and velocities of 
the particles lead to infinitesimal changes in the collision normal and collision 
time. This, in turn, leads to infinitesimal changes in both positions and veloci- 
ties right after the collision. Throughout the paper, primes denote coordinates 
in phase space and tangent phase space just before a collision while non-primed 
quantities refer to coordinates just after the collision. For colliding particles i 
and j, one finds 



(£ + X) 



I -S s 







V 



S I -S 
-Q Q l-S S 
Q -Q S l-S 



5v' 



(6) 



where X and I are the Ad x Ad and dx d identity matrices, respectively, and £ 
is the Ad X Ad collision matrix, which can be written in terms oi dx d matrices 
S and Q specifying the collision dynamics in tangent space (see, for instance, 
Refs. [17]). 

Let Z{t) be the 2dN x 2dN matrix which performs the single-particle trans- 
formations Z{t) for all particles during free flight according to Eq. ([5]). Let 
Lz be the 2dN x 2dN matrix performing the transformations associated with 
collision number z of two particles at time t^, as written in Eq. ([H]), while 
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leaving the other particles unaffected. M^g(t) in Eq. ([T]) then is a product of 
these matrices for the sequence of collisions 1,2, ... ,z occurring between time 
ti = and t and the free flights in between the collisions. Its specific form 
reads 



(t) = Z{t - Q ■ U ■ Z{t, - t,_i) ■ U_i ■ . . . ■ Z{t2 - h) ■ Li . (7) 

In principle, M^,, (t) in Eq. ([1]) can now be derived for finite times for any finite 
system with given initial conditions, after which the eigenvalues of M^g{t) can 
be determined. 



3 Goldstone modes and the extended Enskog equation 



At this point, in order to be able to compare the random-matrix approach 
described in this paper to our previous analytical results obtained by using 
an extended Enskog equation [7], it is useful to review the derivation of the 
extended Enskog equation and the approximations made therein. 

In Ref. [7], it was shown that the sinusoidal modes found in the simulations 
may be interpreted as being Goldstone modes. These occur in systems with a 
continuous symmetry, such as the symmetries associated with the zero modes. 
Translation invariance, for instance, causes the evolution operator to commute 
with the translation operator. As a result, the two operators have a set of 
common eigenfunctions, which are sinusoidal perturbations of the continuous 
symmetry that grow or shrink slowly with time. For these modes to stand out 
among the continuum of other Lyapunov exponents, their wavelength must 
be large compared to the typical length scale of the mean free path. That is. 



vk<^u , (8) 

where v is the average absolute value of the velocity of the particles, and v is 
the average single-particle collision frequency. The modes found in simulations 
are linear combinations of the Goldstone modes for wave vectors of equal 
size belonging to left- and right-moving waves. The qualitative behaviour 
seen in the simulations [12] can be understood entirely from this. Also, the 
inconsistencies in the propagation velocity and amplitude noted in Ref. [19] 
can be completely understood from the behaviour of linear combinations of 
left- and right-moving Goldstone modes. 

The mode corresponding to a particular Lyapunov exponent depends, in gen- 
eral, on the initial conditions of all particles in a way that is far too com- 
plicated to allow exact specification. To find the exponents from analytical 
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calculations, therefore, one has to resort to a statistical approximation. In our 
earlier paper [7] , this was accomplished by making assumptions similar to the 
StoBzahlansatz in the Boltzmann equation. The system is described not as a 
large number of separate particles, but rather by a distribution function of 
position, velocity, etc, which can be used to obtain the average behaviour. 
The pre-coUision pair distribution functions are approximated by a product of 
independent one-particle distributions, which is valid at low densities. In the 
Boltzmann equation, in addition to the StoBzahlansatz, the two one-particle 
distributions are evaluated at the same position r. The Enskog equation is a 
heuristic generalisation of this, in which the pair-distribution is approximated 
by the product of two one-particle distribution functions evaluated at the 
actual positions of the two particles, multiplied by a factor Xe(^) equal to 
the equilibrium pair correlation function at contact evaluated as a function 
of the density n at that point. More details on the Boltzmann equation and 
the Enskog theory of dense gases relevant to this study may be found in 
Refs. [201I2T] . 

To describe the dynamics in tangent space, we have previously derived |7] 
a generalised Enskog equation for the single-particle distribution function 
/(r, V, (5r, 5v, t) for the coordinates of a particle (r,v) and the tangent space 
vectors 5y and 5v. The latter are described by a single-particle distribution 
function which, in the thermodynamic limit, depends smoothly on position, 
velocity, and time, just like the velocity distribution in the ordinary Enskog 
equation. If in addition the distribution functions of the tangent space vectors 
of two particles about to collide are assumed to factorise in a way similar to the 
distribution of their velocities, one ends up with a generalised Enskog equation 
that includes the tangent space variables. If the tangent space variables 5y and 
5v are integrated over, this equation reduces to the standard Enskog equation. 

Because 5r and 5v are infinitesimal, the dynamics in tangent space are linear 
in these quantities. Therefore, from the extended Enskog equation, one may 
obtain closed linear equations for the time evolution of the average first mo- 
ments 5r(r,v,t) and 5v(r,v,t), which are the averages of the tangent space 
vectors 5rj and 5vj of the particles in a small region around position r and 
velocity v at time t. The result is a set of linear equations for the averages [7], 
reading 



— 5v(r, V, t) = -V ■ Tr<^v(r, v, t) + Cs(5v(r, v, t) + Cq5v{v, v, t) . (10) 



The linear functional collision operators Cg and Cq are associated with the 
matrices S and Q, and specified by the collision integrals for a particle with 
outgoing velocity v over all outgoing velocities u for the other particle. 




(9) 
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Cs (5q(r, v,t) 



a--(v-u)<o 



d-l 



CT • (V - u)|0m(u) 



X {(5q(r,v', 



t) + S ■ [ (5q(r + a(T, u', t) - 6q{ 



r,v',t)] 



5q(r, V 



(11) 



CQ(5r(r, v,)!:) 



L 



<T-(V-U)<0 



d-l 



CT ■ (V - u)|0m(u) 



xQ ■ [(5r(r + a(T,u',t) - (5r(r, v',t)] , 



(12) 



where Sq{r, v, t) denotes either 5r(r, v, t) or 5v(r, v, t). As in Sec. [2]and Fig. [H 
the primes denote the variables before the colhsion with colhsion normal <t 
and outgoing velocities v for one particle and u for the other, while 0m (u) is 
the Maxwell distribution. 

The equations for the first moments can be solved by using spatial Fourier 
transforms. The solutions have the form 



where k is the wave vector. For the Lyapunov exponents there is reasonable 
quantitative correspondence to the results of simulations [7]. Surprisingly, 
it turned out that the Stofizahlansatz still affects the leading order in the 
density [71122] . 

The extended Enskog equation, the equation for the first moments, and its 
solutions are discussed in more detail in Ref. [7]. 



4 Products of random matrices and the weak disorder expansion 

Because Eq. (JTj) contains a product of similar matrices, a more natural ap- 
proach to the problem of Lyapunov exponents in many-particle systems may 
be to consider the tangent-space maps as products of random matrices. Some 
assumptions regarding the randomness of the generating dynamics must be 
made to replace the matrices by a random ensemble, similarly to the assump- 



tions used in, for example, Refs. pi23f24f25] ). The matrices in the ensemble 




must be symplectic, because of the time-reversal symmetry of the dynamics, 
yet satisfy quite strong additional restrictions and correlations between their 
elements. In the case of a general many-particle system, it is not directly 
obvious how to deal with these remaining correlations, as pointed out in 



Attempts to describe the Goldstone modes through random-matrix theory 




(13) 




Ref. |6]. 
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have been made by Eckmann and Gat [15], and by Taniguchi, Dettmann, 
and Morriss [T6|[T7]. These works considered quasi-one- dimensional systems, 
in which particles cannot move past each other, and any particle collides only 
with its two neighbours. As a result, the calculations are very much simplified, 
and indeed fail to produce quantitatively reliable results for fully two- and 
three-dimensional systems. It should be noted in this context that for the 
same reason quasi-one-dimensional systems do not satisfy the StoBzahlansatz, 
which does not hold if pairs of particles that have already collided collide again. 
Additionally, these calculations make use of approximations that are quite 
similar to the weak-disorder expansion, while the applicability of the weak- 
disorder expansion is uncertain. In the present section, I will indicate how the 
requirements of the weak-disorder expansion correspond to the assumptions 
needed to derive the extended Enskog equation and the equations for the first 
moments in the tangent space of a single particle, Eqs. and (ITU]) . 

For the problem under consideration the matrix product M^g()f:) is given in 
Eq. ([7]). The matrices can be taken together in pairs corresponding to a 
collision j, Lj, and following free flight, Z{tj+i — tj), to produce a product 
of correlated sparse random matrices 



After a collision, following a velocity-dependent free-flight time, each particle 
will be involved in another collision, with the outgoing velocity of the previous 
collision as the new incoming velocity. The matrices describing such collisions 
are correlated. 

In the next two paragraphs, the weak-disorder expansion is introduced, along 
with the conditions that must be met. In the remainder of the section, it is 
shown how these requirements can be met for the product in Eq. (|T4|) and 
how the Lyapunov exponents can be derived. 

In preparation for this, let us consider a quite general product of z matrices 
of the form 



in which e is a small perturbation parameter and A is the average matrix, 
assumed to be diagonalisable. The matrix M thus expands to 




(14) 
(15) 



M = A, ■ A,_i 
Aj = A + eXj , 



• . . . ■ A 



(16) 
(17) 



A'-^ ■ X, ■ A-'- 



+ . . . 
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The terms of higher orders in e are not necessarily smaller than the first-order 
terms, because the number of terms grows faster than the inverse of their size 
for z — )■ oo. However, if the are independent of one another this problem 
disappears and the sum converges [2^1177] . 

Let Xj be independent and identically distributed, while every element has 
zero average. Derrida, Mercheri, and Pichard |28j have derived that the growth 
rates Aj of the eigenvalues of the product can to first order in e be expressed 
in the eigenvalues Ki of the average matrix A, provided these eigenvalues exist 
and are non-degenerate. More specifically. 



Here, is the eigenvector of A belonging to the eigenvalue Ki. This theorem 
was later extended to the degenerate case by Zanon and Derrida [27]. In the 
following sections, it will be shown that the eigenvalues of A in the present 
problem can be chosen arbitrarily, to be distinct, but approaching unity, so 
as to simplify the derivation. For more details on random matrices and the 
weak-disorder expansion, see Crisanti, Paladin, and Vulpiani [26j. 

In order to apply the weak-disorder expansion to the problem under consider- 
ation, the matrix product M^g(t) in Eq. (IT^ must first be rewritten into the 
form of Eq. (fT6|l . i.e., in such a way that the matrices become independent 
and that an expansion parameter e as introduced in Eq. (fTTI) exists. 

Instead of calculating the spectrum of Lyapunov exponents of the original 
product M^Q(t), the spectrum of a very similar product will therefore be 
considered, namely 



■ M^,{t) ■ Oo = ■ B, ■ 0,_i ■ 0;i, ■ B,_i ■ . . . ■ Bi ■ Oo , (20) 



where Oj may depend on the entire phase space and may be chosen in any 
suitable way. Without loss of generality, Oq, may be taken to be unitary, and 
the other Oj are chosen in such a way that the 2dN x 2dN matrices Aj, defined 




+ . . . 



(19) 



by 
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(21) 



satisfy the requirements of the weak-disorder expansion. The matrices Oj are 
transformations to another orthogonal basis just before and back after every 
colhsion and free flight, while the matrices Aj are the tangent space maps in 
these new bases. The basis vectors make up the columns of the matrix Oj. 

In Sec. [5], I derive an equation for the components of these vectors, and 
consequently for the elements of Oj, from the requirements of the weak- 
disorder expansion. These basis vectors may depend on the position of the 
system in phase space, i.e., the positions and velocities of the particles. The 
matrices Oj need not be unitary. Any kind of invertible transformation will 
do, as long as the transformation of a product of tangent space maps is equal 
to the product of the transformed maps. Here, however, I restrict myself to 
linear transformations, because they suffice for the present purpose. 

The Lyapunov spectrum of the product 0^^ ■ M^y(t) ■ Oq, Eq. fEUl) . can then 
be related to the Lyapunov spectrum of the original product M^o(t). As 0^ is 
a transformation from one orthogonal basis to another, it can be written as a 
product of a diagonal matrix Dj and a unitary matrix Uj, 



The first diagonal matrix Do is the 2dN x 2dN identity matrix. The Lyapunov 
spectrum of D-^ ■ \J~^ ■ M^^ (t) ■ Oq can be determined using the Osledec theorem, 
which states that the growth of an arbitrary i dimensional volume is dominated 
by the i largest Lyapunov exponents. The unitary transformations and 
Oo do not affect the volume mapped by the product, and so the growth of a 
volume mapped by LJ-^ ■ M^g(t) ■ Oq is the same as that of a volume mapped 
by M^g(t). Let Ui be the i-th largest eigenvalue of D-^. One can write D^^ 
as a product of z diagonal matrices and apply the Osledec theorem to this 
product as well. The Lyapunov exponents of this product are equal to \ia{ui)/z. 
By starting from the volume mapped by IJ-^ ■ M^^it) ■ Oq and applying the 
Osledec theorem to its map by D^, one obtains that the growth rate of an i 
dimensional volume mapped by the entire product in Eq. (!20l) . is equal to the 
sum of the largest i Lyapunov exponents of each of the two products M^g(t) 
and D-i. The Lyapunov exponents Aj of the complete product in Eq. (120]) can 
therefore be written as 



0, = U, ■ D, . 



(22) 



— Aj + 



In ooi 



(23) 



z 
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5 Basis vectors for the transformations to independent matrices 

In the previous section, a framework was set up which allows one to implement 
the weak-disorder expansion. In order to materialise this, time- dependent 
matrices Oj must be constructed that make the random matrices Aj, Eq. (121 p . 
occurring in the product Eq. fl20|) independent. The columns of Oj are made up 
of a time-dependent set of 2 dA^- dimensional basis vectors Ox{t), enumerated 
by X which runs from 1 to 2dN. It should be noted that it is convenient to 
choose vectors Ox{t) which are orthogonal. 

In the many-particle system under consideration, the correlation between the 
matrices derives from the fact that information is carried between collisions 
through the position and velocity of a particle. The basis vectors in which 
to express the tangent space maps Aj must be chosen in such a way that 
this correlation is removed. In this section, this is achieved by focussing on 
the subsequent collisions of one particle. Most correlation with collisions of 
other particles can be removed by using the StoBzahlansatz. The remaining 
correlation, which involves only a small number of particles can be used to 
determine a set of equations for Oxit). 

Let Pi^p be the product of matrix Aj_i ■ . . . ■ Aj+i • Aj which describes the 
dynamics in tangent space of all the particles at collision i involving particle 
p and during the subsequent free flight of particle p until it collides with 
particle q in collision j (see Fig. [2]). Using Eq. fl?T]) . the x,y element of Pjp 
can be written in terms of the basis vectors Ox{t~) and Oy{t~) as 

(P.p).. = 1^1^ ■ B,-i . . . B, ■ Oy{t-) , (24) 

where the superscript in t~ etc. denotes a time an inflnitesimally small interval 
before the collision at that time. Because of the StoBzahlansatz, any collisions 
of particles with no common history with p or g are not correlated with 
collisions of p or q. As our purpose here is to remove all correlation, and 
these collisions are already uncorrected, they do not need to be considered 
explicitely. 

In order to remove the remaining correlation, and for the matrices kj to 
become fully independent and to satisfy the requirements of the weak-disorder 
expansion, one must consider the remaining correlation, which is due to cor- 
relations between particles which share a common history. Though there will 
be correlation between collisions of such particles with other particles, the 
corresponding tangent space maps Aj should be independent. For this, it is 
suflicient that the product P^^p be independent of both ^ and Pj^q. This 
ensures that the matrix products are independent of the history of a particle 
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position 




Fig. 2. Schematic representation of the collisions and free flights of particle p and the 
particles it interacts with at collisions i, j, k, and, indirectly, I, at times U, tj, tk, and 
ti, along with the corresponding tangent space maps and their products. Because 
of the Stofizahlansatz, the only remaining correlation between random matrices in 
the product is due to the information carried by particles between collisions. If the 
matrices Aj, . . . , A;_i are constructed in such a way as to make Pj^p independent 
of Pj^p and Pj^q, then this correlation is eliminated. 



before its most recent collision, and consequently that collisions of any other 
particles which share history with p and q can be ignored. The component 
of Ox{t) belonging to a specific particle should not depend on anything but 
its history since its most recent collision. Just before collision j, therefore, 
only the recent history of particle p can lead to correlation. Consequently, 
only contributions to {Pi,p)xy from the projections onto the tangent space of 
particle p can lead to correlations with P^-^p and Pj^g, i. e., we may restrict 
ourselves to the term 



■^^).p,.B,.....B..o.(,r). (25) 



where Vp denotes the projection of the 2 diV- dimensional tangent space onto 
the 2d-dimensional tangent space of particle p. Furthermore, we may write 
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. . . ,5rj^^(rAr, VAr,t),5vj^^(rAr, VAr,t)) , (26) 

where 5r^^'*(r, v, t) . . . 5vj^^(r, v,t) are (i-dimensional vector functions which 
are different for every x. As the particles are identical, the functions 5r^^^(r, v, t), 
6y2'\v, V, t), etc are the same for all particles, as are 6\i'\y, v, t), 6\2^\y, v, t), 
etc., so that the indices can be dropped and they can simply be writen as 
5r*^^^(r, V, t) and 5v'^^^(r, v, t). 

As the only information carried by particle p between two collisions is its 
position Yp and velocity Vp, the product Pi^p is independent of Pj^p and Pj^q if 
Pi^p is independent of the position and the velocity of the particle p between the 
two collisions i and j. To this end, the functions 5y^^\y^ v, t) and 5\^^\y^ v, t) 
must be chosen in such a way that Pj^p becomes independent of and Vp. 
For this it is necessary that at least the average (Pi,p) is independent of the 
velocity of particle p after collision where the average (.) is calculated over all 
parameters, including the velocity of all particles that particle p has collided 
with, but excluding the position and velocity of particle p after collision i. If 
this condition is met, averages of products of first moments of the matrices 
Pj^p with those of Pj p and Pj^q behave as if the matrices are independent. 



(Pj,pPj,p) — (Pj,p) (Pj,p) , (27) 
(P.,pP,,,) = (P.,p)(P,g . (28) 

The condition of independent first moments is a necessary prerequisite for the 
independence requirement of the weak-disorder expansion, if not necessarily 
sufficient. It is assumed here that the remaining correlation can be removed 
by further refining the maps Oj [see Eq. (12T!) ]. if necessary through the use 
of nonlinear maps. As only the higher moments of the matrices would be 
involved in this, it will not produce extra contributions to the leading order 
in the perturbation parameter in Eq. (fT9l) . 

One must consider the relevant contribution to (Pj,p) from particle p in Eq. ( !25l) . 
which is proportional to 



((5r(^')(rp,vp,t-),5v(^)(rp,vp,t-)) ■Pp-B,_i...B,-o,(tr)) , (29) 

where Vp denotes the projection of the 2 dA?"- dimensional tangent space onto 
the 2(i-dimensional tangent space of particle p. The functions 5r*^^^(rp, Vp, t) 
and 5v'^^)(rp, Vp, t) do not contain any of the quantities which are averaged 
over. Therefore the quantity 
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(5r(^)(rp, v„ tj), 5v(^)(r„ v„ tj)) ■ {V, ■ B,_i . . . B, ■ o,(tr)) 



(30) 



must be independent of Vp and Vp. If the average has the same functional de- 
pendence (up to an arbitrary prefactor) on the position and velocity of particle 
p just before collision i as dv^^^Vp, Vp, tj), then the innerproduct in Eq. fl5U]) is 
equal to an arbitrary constant times the norm of (Sr^^\rp, Vp, t), 5\'^^\rp, Vp, t)) 
squared. In this case, Pj^p represented in the two bases sets just before collisions 
j and i is independent of the position and velocity of particle p between the 
two collisions. 

Consequently, one arrives at an eigenvalue-type equation for the tangent-space 
basis vectors with a strictly positive, but otherwise arbitrary, scalar prefactor 
V, which may be different for each combination of basis vectors. The collision 
parameters must be averaged out, and their distribution is the same as in 
Enskog theory. Let 5r(r,v,t) and 5v(r,v,t) be shorthand for (5r^^^(rp, Vp, t) 
and 5\'^^\vp,\'p,t), the eigenvalue equation reads 



where the collision normal is &, and the incoming velocities v' and u' of two 
particles are mapped to the outgoing velocities v and u, respectively. The 
4d X 4d matrix Cj represents the collision dynamics at collision j as given in 
Eq. ([6]) and the 2d x 2d matrix Z{t) defined in Eq. ([5]) describes the free- 
flight dynamics. The time between the two collisions of the particle is denoted 
by r. The velocity-dependent average collision frequency is represented by 
z/(v) = l/f(v), where r(v) is the average mean free time as a function of 
the particle velocity. The various solutions of Eq. (13T!) correspond to different 
basis vectors Ox{t), which should be orthogonal. 

The strictly positive time-dependent prefactor V controls the eigenvalues of 
A as well as the growth rate of the eigenvalues of Dj in Eq. (!22|) . In principle, 
it can be chosen in any arbitrary way, as long as it is strictly positive, as any 
growth can be accounted for either in the basis vectors or A. The choice of 
V does not affect the final result for the Lyapunov exponents, as can be seen 




5r{r — TV, v', t — t) 



xZ{t) ■ 




(31) 
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from Eq. ( 123|) . If V is taken to be unity, A is equal to the identity matrix, and 
the eigenvalues of A are all equal to unity. All growth will be accounted for 
in the time-dependence of the functions (5r(r,v,t) and (5v(r,v,t). The weak- 
disorder expansion, however, has only been proven for non-degenerate systems 
and systems with doubly degenerate eigenvalues [27j. I therefore choose V 
differently for every basis vector, so that the eigenvalues of A become distinct 
and the weak-disorder expansion holds. Furthermore, P is chosen so close to 
unity that terms proportional to P — 1 become negligible compared to other 
terms. This greatly simplifies the calculation. 

Let T be the operator /q°° druiy) exp[— rz/(v)]. By substituting Eqs. (ITTi) and 
(fT2l) into Eq. ([31]), one finds 



-5r(r, V, t) = T {[r(v)Cs + 1 + rr(v)CQ] (5r(v, r - rv, t - r) 

+ r[f(v)Cs + l]5v(v,r-rv,t-r)} , (32) 

— 5v(r, V, t) = T {f (v)CQ5r(v, r - rv, t - r) 

+ [f(v)Cs + l](5v(v,r-rv,t-r)} . (33) 

The operator T equals the unity operator when working on functions that do 
not depend on r, and therefore one may multiply the left side of the equations 
by it. As both sides then contain the same integral over r, a change of variables 
may be performed, namely r' = r — rv, t' = t — r. After renaming the variables 
the equations become 



— T5r(r + rv, v, t + r) = [r(v)Cs + l]5r(v, r, t) 

+ Trf (v) [CQ^r (v, r, t) + (Cs + l)5v(v, r, t)] , (34) 
^T6^f{r + tw, w,t + T)= r(v)CQ(5r(v, r, t) + [r(v)Cs + l]5v(v, r, t) .(35) 



In the continuum limit the mean free time and the mean free path become 
small in comparison with the typical scales of the solutions. The quantities 
5y{y — rv, V, t — r) and 5v(r — rv, v, t — r) can thus be expanded around 
5r(r, v,i(:) and 5r(r, v,t), and all terms of higher order in r can be neglected. 
After substitution of Eq. (135|) into Eq. and performing the integration 
over r, one finds 
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= [r(v)Cs + l]5r(v, r, t) + f (v)5v(v, r, t) , 
= f (v)CQ5r(v, r, t) + [f (v)Cs + l](5v(v, r, t) . 



(36) 



(37) 



Rearranging the terms and dividing the equations by r(v), followed by taking 
the limit of P — ?■ 1, leads to two equations very similar to Eqs. and f lTU]) . 



d d 

— 5r(r, V, t) = -V ■ ^5r(r, v, t) + (5v(r, v, t) + Cs(5r(r, v, t) 
OT ar 



+ {V-1) 



Cs + 



r(w)^ 



5r(r, V, t) + 5v(r, v, t) 



d d 

— (5v(r, V, t) = -V ■ ^(5v(r, v, t) + Cs(5v(r, v, t) + CQ5r(r, v, t) 
at or 



+ {V-1] 



CQ(5r(r,v,t) + Cg + 



(5v(r, v,t) 



(38) 



.(39) 



As V has been chosen arbitrarily close to unity, and the other terms are 
nonzero, the terms containing V — 1 are negligible. We are then left with equa- 
tions identical to Eqs. and f llOp . The latter were derived from the Enskog 
equation by using the same approximations in a different order. The equations 
can be Fourier transformed and have solutions of the form of Eq. f|T3|) . The 
solutions give the components of the basis vectors that make up Oj in Eq. fl2T]) . 
As the average of the transformed matrices. A, is also arbitrarily close to unity, 
to leading order the Lyapunov exponents are simply equal to the growth rates 
of the basis vectors, the solutions of Eqs. ([9]) and ( ITOi) . 



6 Perturbation parameter 



The pursuit for independently distributed matrices Aj has yielded the same 
equations as were found in Ref. [7j. However, in order for A^ to fully satisfy 
the requirements of the weak-disorder expansion, and thus for the Lyapunov 
exponents to equal the A found from Eq. ([9]) and ( fTOl) . Aj must be of the form 
of Eq. (ITTl) . A small perturbation parameter e must exist, i.e., the higher-order 
corrections in Eq. (fT9|l must be small compared to the leading order. 

The solutions to Eqs. ([9]) and (ITOl) . which are of the form given by Eq. (fT3|) . 
describe the vectors that make up the columns of the matrix 0^ associated 
with a specific collision b at time U and the subsequent free flight.. We may 
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estimate the elements e)C^ of the corresponding disorder matrix by considering 
the inner product of these basis vectors given by Eq. (fT3|) . with the tangent- 
space colhsion-matrix and free-flight matrix in between. Here i and j determine 
not only the wave vectors, but also the zero mode modulated by them. The 
contribution from particle I to the tangent space vector at time can be 
written as [7] 

5q«(r^ h) = [Aq«(vO] exp(ik, • + X^b) , (40) 
Aq«(vO=Aoq«(vO + of^) , (41) 



where Aoq*-*^ denotes a linear combination of the zero modes that belongs 
to the mode of Aj. Note that the components of AQ^q(v;) are of order 1/A^ 
due to normalisation. If fcj becomes large, then, because of normalisation, the 
components are of order 1/A^ < ki/N. Note that modes with different wave 
vectors kj decouple on average after taking the continuum limit. 

For all particles in the time interval Atb ~ 1/{jjN) between this collision b, 
and the next one of the entire system, 6+1, the position changes by a small 
amount vjAtfo. More explicitly, the matrix A;, related to any specific collision 
b and the directly following free flight [see Eq. ( l20l) ] can be expressed as 

... E;iiEq=r,v lAh+i-h) ■ U-(5qW(r,,v,,t,)] •^q(^-)(r;,v;,t,) 
^ '''' E^IiE,=r.v5q(^Kr;,v,',t,)-5qO)(r;,v;,t,) 

X exp[(-iki ■ v; + Xi)Ah] , (42) 



Out of the total of particles, only two nearby particles are involved in 
each collision. To leading order in vki/{yN), their perturbations are the zero 
modes, which map onto themselves under collisions. All the other particles are 
unaffected by the collision matrix. As a result, after the short time-evolution 
after collision b the modes are mapped exactly onto themselves, except for 
small terms for each particle due to the free flight, and larger collision terms 
for the two colliding particles. 

The expressions in Eqs. (HOj) and fHT]) may be inserted into Eq. (142|) to 
estimate the various terms. Upon inserting the orthogonality of the basis and 
subtracting A, which approaches unity, one finds 



One may substitute this result into Eq. (|T9l) while considering that the eigen- 
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values Ki of the average tangent space evolution matrix A are of order one. In 
this way, it is found that the higher-order terms in the expansion are all of 
order [vkiY / {p"^ N) . 

In other words, the slowly growing and decaying modes couple sufficiently 
weakly to the other modes if vk/v ^ 1. This is the case if the modes vary 
slowly in space, i.e., if the wavelength is smaller than the mean free path, as 
in Eq. ([8]). In this perturbation expansion all the higher-order corrections in 
Eq. (I19p are of order {vk/vY, one order oivk/v smaller than the leading-order 
term. The perturbation parameter needed for the weak-disorder expansion is 
therefore vk/v, and not vk/{uN), because of the size of the disorder matrices. 
As the wavenumber k appears in the matrix 0^, it appears in A;,, and can 
appear in the perturbation parameter. 

Note that the perturbation expansion can only be used to calculate the small 
Lyapunov exponents, when vk/u < 1. If k becomes too large, the remaining 
terms in Eq. f|T9l) may not converge. 



7 Comparison and discussion 

Exactly the same approximations are used in both the Enskog approach 
(Sec. [3]) and the random-matrix approach (Sees. IH and [H]), namely the 
continuum limit, the StoBzahlansatz, and the approximation of long wave- 
length, Eq. ([8]). The same results are found. However, the method for arriving 
at these results is quite different, and the two approaches shed different light 
on the approximations and their validity. 

In the derivation of the extended Enskog equation using kinetic theory, one 
starts with the approximations made in the standard Boltzmann and Enskog 
equations, namely the StoBzahlansatz and the continuum limit, keeping in 
mind that the wavelength of the modes must be longer than the mean free 
path. Next, the equations are multiplied by the first moments and integrated. 

In the random-matrix approach, only a very general assumption of random- 
ness is made at the beginning. Then weak disorder is used, along with the 
approximation of long wavelength, followed by the Stofizahlansatz, which is 
needed to simplify the equations for the basis. And only then the continuum 
is limit taken. 

As to the Stofizahlansatz, in the Enskog approach it is used earlier, and the 
consequences of the Stofizahlansatz are clearer. In Ref. [7j we showed that 
this approximation may not be sufficient for finding the Lyapunov exponents 
to leading order in the density. From the kinetic-theory approach it is clear 
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how to add correction terms to the equations for ring coUisions and other 
similar trajectories. In the random-matrix approach this is not so transparent, 
as the weak-disorder expansion theorem simply requires independence of the 
matrices. It is not directly obvious how to extend the theorem to products of 
correlated matrices, or how to include ring-collisions in the equation for the 
basis sets for which the matrices become independent. An exploration of the 
ring-collision terms within the framework of the kinetic-theory approach was 
made in Ref. [22] . 

In the derivation of the generalised Enskog equation it is assumed, because 
the tangent space is linear, that the behaviour of the system can be described 
by the average behaviour of the first moments of the tangent space vectors. 
This is equivalent to the requirement that the averages of the first moments of 
the matrix elements behave as if the matrices are independent. Although an 
expansion in vk/v of the type written down in Eq. f[T^ was not explicitly 
considered in the kinetic-theory approach, the same requirements of long 
wavelength and the StoBzahlansatz were found to be needed. Clearly, the 
original assumption made when deriving the generalised Enskog equation is 
not as trivial as it may seem. 

The continuum limit is taken at the very end in the random-matrix approach. 
Because of this, it is much more transparent than in the generalised Enskog 
equation how to add terms for finite-size effects to Eqs. ([9]) and (fTOj) . In par- 
ticular, the longitudinal and transverse modes, and also modes with different 
values of k, are not perfectly orthogonal in systems with a finite number 
of particles. Such nonorthogonality is directly related to fluctuations (in for 
instance the energy density) and its decay is related to the decay of fluctuations 
and correlation of position, energy, and momentum. The average fluctuations 
of the inner products between modes with different wave numbers are zero, 
and therefore finite-size effects cannot contribute to the Lyapunov exponents 
to the leading order in k. However, the relative size of these fluctuations is 
of order l/\/N and the total number of wave vectors is of order A^. There 
could be contributions to the higher-order terms on the right side of Eq. f[T^ . 
which are of order k'^/N, leading to corrections to the Lyapunov exponents of 
order k"^. 



8 Conclusions 

The Lyapunov exponents close to zero can be related to Goldstone modes. In 
this paper, I have shown that the equations for the small positive Lyapunov ex- 
ponents found in Ref. [7] by using Enskog theory can be derived from random- 
matrices by use of the weak-disorder expansion. Precisely the same approx- 
imations are necessary in both approaches, namely the StoBzahlansatz, the 
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continuum limit, and the approximation of long wavelength. In the random- 
matrix approach these approximations are needed in order to find a suit- 
able description of the basis to meet the requirements of the weak-disorder- 
expansion theorem. The approaches were thus found to be equivalent. 

The different order in which the approximations are made, however, enables 
us to gain understanding into the consequences of the approximations and 
the uncertainties introduced by them. In the kinetic-theory approach, it is 
more transparent how to add corrections for ring collisions and correlations, 
whereas in the random-matrix approach it is more clear how to add corrections 
for shorter wavelengths and finite-size effects. The latter may give insight into 
connections between the Lyapunov exponents and the decay of correlations 
and fluctuations. Similar connections were suggested by Taniguchi and Mor- 
riss [29]. Further investigation of the continuum limit and the validity of the 
weak-disorder expansion could perhaps also provide insight into the precise 
nature of the destruction of the Goldstone modes as the wavelength becomes 
too short compared to the mean free path. 

In systems with soft potentials, the calculations described in this paper and in 
Ref. [7j become more complicated. Though the Goldstone modes are present 
in simulations of gases of soft particles, the degeneracy disappears [TT|18|30] . 
This could be related to the fact that the zero modes are not as trivial in 
soft-potential systems as they are in system with hard-core interaction |22]- 
Soft potentials can produce qualitative differences in the Lyapunov exponents, 
as has also been concluded in Refs. [3T1I321. 
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